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Factorization Breaking in Diffractive Photoproduction of Dijets 
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Abstract. We have calculated the diffractive dijet cross section in low-Q'^ ep scattering in the HERA 
regime. The results of the calculation in LO and NLO are compared to recent experimental data of the 
HI collaboration. We find that in LO the calculated cross sections are in reasonable agreement with the 
experimental results. In NLO, however, some of the cross sections disagree, showing that factorization 
breaking occurs in that order. By suppressing the resolved contribution by a factor of approximately three, 
good agreement with all the data is found. The size of the factorization breaking effects in diffractive dijet 
photoproduction agrees well with absorptive model predictions. 
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1 Introduction 

Diffractive 7p interactions are characterized by an out- 
going proton of high longitudinal momentum and/or a 
large rapidity gap, defined as a region of pseudo-rapidity, 
rj = — lntan0/2, devoid of particles. It is assumed that the 
large rapidity gap is due to the exchange of a pomeron, 
which carries the internal quantum numbers of the vac- 
uum. Diffractive events that contain a hard scattering are 
referred to as hard diffraction. A necessary condition for 
a hard scattering is the occurrence of a hard scale, which 
may be the large momentum transfer in inclusive deep- 
inelastic ep scattering, the high transverse momentum of 
jets or single hadrons, or the mass of heavy quarks or of 
W-bosons produced in high-energy 7p, ep or pp collisions. 

The central problem in hard diffraction is the ques- 
tion of QCD factorization, i.e. the question whether it 
is possible to explain the observed cross sections in hard 
diffractive processes by a convolution of diffractive par- 
ton distribution functions (PDFs) with parton-level cross 
sections. 

The diffractive PDFs have been determined by the 
HI collaboration from a recent high-precision inclusive 
measurement of the diffractive deep inelastic scattering 
(DIS) process ep eXY, where F is a single proton or 
a low mass proton excitation pQ. The diffractive PDFs 
can serve as input for the calculation of any of the other 
diffractive hard scattering reactions mentioned above. For 
diffractive DIS, QCD factorization has been proven by 
Collins 0. This has the consequence that the evolution 
of the diffractive PDFs is predictable in the same way 
as the PDFs of the proton via the DGLAP evolution 
equations. Collins' proof is valid for all lepton-induced 



collisions. These include besides diffractive DIS also the 
diffractive direct photoproduction of jets. The proof fails 
for hadron-induced processes. 

As is well known, the cross section for the photopro- 
duction of jets is the sum of the direct contribution, where 
the photon couples directly to the quarks, and of the re- 
solved contribution, where the photon first resolves into 
partons (quarks or gluons) , which subsequently induce the 
hard scattering to produce the jets in the final state. So, 
the resolved part resembles hadron-induced production of 
jets as for example in pp collisions. Dijet production in 
single-diffractive collisions has been measured recently by 
the CDF collaboration at the Tcvatron 3 . It was found 
that the dijet cross section was suppressed relative to the 
prediction based on older diffractive PDFs from the HI 
collaboration 0] by one order of magnitude |3] . From this 
result we would conclude that the resolved contribution in 
diffractive photoproduction of jets should be reduced by 
a correction factor similar to the one needed in hadron- 
hadron scattering f3| . This suppression factor (sometimes 
also called the rapidity gap survival probability) has been 
calculated using various eikonal models, based on multi- 
pomeron exchanges and s-channel unitarity The di- 
rect and the resolved parts of the cross section contribute 
with varying strength in different kinematic regions. In 
particular, the x-y-distribution is very sensitively depen- 
dent on the way how these two parts of the cross sec- 
tion are superimposed. Near ~ 1 the direct part dom- 
inates, whereas for < 1 the resolved part gives the 
main contribution. However, in this region also contribu- 
tions from next-to-leading order (NLO) corrections of the 
direct cross section occur. Therefore, to decide whether 
the resolved part is suppressed as compared to the exper- 
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imental data, a NLO analysis is actually needed. This is 
the aim of this paper. For our calculations we rely on our 
work on dijet production in the inclusive (sum of diffrac- 
tive and non-diffractive) reaction 7 + p — > jets + ^'i^ [Z], 
in which we have calculated the cross sections for inclu- 
sive one-jet and two-jet production up to NLO for both 
the direct and the resolved contribution. The predictions 
of this and other work ,8 have been tested now by many 
experimental studies of the HI and ZEUS collaborations 
Very good agreement with the experimental data 
pillUj has been found. From these comparisons it follows 
that a leading order (LO) calculation is not sufficient. It 
underestimates the measured cross section by up to 50% 

The question whether the resolved cross section needs 
a suppression factor, can be decided first by looking at 
the shape of those distributions which are particularly 
sensitively dependent on the resolved contributions, as 
for example the a;^-distribution for the smaller x-y or the 
£'T-distributions at small Et- Because of the interplay 
of direct and resolved contributions, LO calculations are 
not sufficient, in particular, since the NLO corrections are 
much more important for the resolved than the direct part. 
This is even more important if one looks at the normal- 
ization of the differential cross sections. 

Recently the HI collaboration |T21 have presented data 
for differential dijet cross sections in the low-|t| diffractive 
photoproduction process ep — *■ eXY, in which the photon 
dissociation system X is separated from a leading low- 
mass baryonic system y by a large rapidity gap. Using the 
same kinematic constraint as in these measurements we 
shall calculate the same cross section as in the HI analysis 
up to NLO. By comparing to the data we shall try to find 
out, whether or not a suppression of the resolved cross 
section is needed in order to find reasonable agreement 
between the data and the theoretical predictions. 

The outline of this work is as follows. In Sec. El we 
specify the kinematic variables used in the analysis and 
describe the input for the calculation of the diffractive di- 
jet cross section. In Sec. 13 we report our results and dis- 
cuss our findings concerning the suppression factor for the 
resolved contributions. Section 0] contains our conclusions 
and the outlook to further work. 



2 Kinematic Variables and Diffractive Parton 
Distributions 

2.1 Kinematic Variables and Constraints 

The diffractive process ep — > eXY, in which the systems 
X and Y are separated by the largest rapidity gap in the 
final state, is sketched in Fig.Q] The system X contains at 
least two jets, and the system Y is supposed to be a proton 
or another low-mass baryonic system. Let k and p denote 
the momenta of the incoming electron (or positron) and 
proton, respectively and q the momentum of the virtual 
photon 7*. Then the usual kinematic variables are 




\ w 



Fig. 1. Diffractive scattering process ep — > eXY, where the 
hadronic systems X and Y are separated by the largest rapidity 
gap in the final state. 



We denote the four-momenta of the systems X and Y by 
px and pY- The HI data ^21 are described in terms of 



Mi and t={p-pY)^ 
r2 qip-py) 



My — Py and xp 



qp 



(2) 



where Mx and My are the invariant masses of the systems 
X and Y, t is the squared four-momentum transfer of 
the incoming proton and the system Y, and Xp is the 
momentum fraction of the proton beam transferred to the 
system X. 

The exchange between the systems X and Y is sup- 
posed to be the pomeron IP or any other Regge pole, 
which couples to the proton and the system Y with four- 
momentum p — Py- The pomeron is resolved into partons 
(quarks or gluons) with four- momentum v. In the same 
way the virtual photon can resolve into partons with four- 
momentum u, which is equal to q for the direct process. 
With these two momenta u and v we define 



pu 

— and zp 
pq 



qv 



qip-py)' 



(3) 



and y 



qp 
kp 



(1) 



x~f is the longitudinal momentum fraction carried by the 
partons coming from the photon, and zp is the corre- 
sponding quantity carried by the partons of the pomeron 
etc., i.e. the diffractive exchange. For the direct process we 
have x^ = 1. The final state, produced by the ingoing mo- 
menta u and V, has the invariant mass Mi 2 = y/ {u + w)^, 
which is equal to the invariant dijet mass in the case 
that no more than two hard jets are produced, q — u and 
P — Py — V are the four-momenta of the remnant jets pro- 
duced at the photon and pomeron side. The regions of the 
kinematic variables, in which the cross section has been 
measured by the HI collaboration ^^li are given in Tab.^ 
With the same constraints we have evaluated the theoret- 
ical cross sections. 

The upper limit of xp is kept small in order for the 
pomeron exchange to be dominant. In the experimental 
analysis as well as in the NLO calculations, jets are defined 
with the inclusive fcT-cluster algorithm with a distance 



M. Klasen, G. Kramer: Factorization Breaking in Diflractive Photoproduction of Dijets 



3 



Table 1. Regions of kinematic variables. 



0.3 < 


y 


< 


0.65 






< 


0.01 GeV^ 






> 


5 GeV 






> 


4 GeV 


-1 < 


'/lab 


< 


2 




XIP 


< 


0.03 




My 


< 


1.6 GeV 






< 


1 GeV^ 



parameter d = 1 |^ in the laboratory frame. At least two 
jets are required with transverse energies Ei^*^ > 5 GeV 
and £;j?*^ > 4 GeV. They are the leading and subleading 
jets with —1 < rjll^'^ < 2. The lower hmits of the jet 
Et^s are asymmetric in order to avoid infrared sensitivity 
in the computation of the NLO cross sections, which are 
integrated over Et |14| . 

In the experimental analysis the variable y is deduced 
from the energy E'^ of the scattered electron y = 1 — 
E'^/Eg. Furthermore, sy = = [q + p)'^ = {px + Py)^- 
xp is reconstructed according to 



2En 



(4) 



where Ep is the proton beam energy and the sum runs over 
all particles (jets) in the X-system. The variables M12, x.^, 
and zp are determined only from the kinematic variables 
of the two hard leading jets with four-momenta and 
pi''*^. So, 

where additional jets are not taken into account. In the 
same way 



Jots 



Ejcts(£^-p.) 



2yE, 



and 



Jots _ 

^-p ~ 



IxpEp 



(6) 



The sum over jets runs only over the variables of the two 
leading jets. These definitions for x^ and zp are not the 
same as the definitions given earlier, where also the rem- 
nant jets and any additional hard jets are taken into ac- 
count in the final state. In the same way Mx can be es- 
timated by M|- = Ml^/i^z^p^'x^^^''). The dijet system is 
characterized by the transverse energies i?^*^ and i?;^*^ 
and the rapidities in the laboratory system yyj^" and 'q\'^^ . 
The differential cross sections are measured and calculated 
as functions of the transverse energy E}^^^ of the leading 
jet, the average rapidity rf°^^ = (?7j^" + '?iab^)/2, and the 
jet separation |Z\7^-'''*'*| = l^yja" — ''f\ah\^ which is related to 
the scattering angle in the center-of-mass system of the 
two hard jets. 



2.2 Diffractive Parton Distributions 

The diffractive PDFs are obtained from an analysis of 
the diffractive process ep — s- eXY , which is illustrated in 



Fig.Q where now is large and the state X consists of all 
possible final states, which are summed. The cross section 
for this diffractive DIS process depends in general on five 
independent variables (azimuthal angle dependence ne- 
glected): Q"^ , X {or (3), Xp, My, and t. These variables are 
defined as before, and x = Q^/{2pq) = Q^/iQ^ + W^) = 
xp(3. The system Y is not measured, and the results are 
integrated over —t<l GeV^ and My < 1.6 GeV as in 
the photoproduction case. The measured cross section is 
expressed in terms of a reduced diffractive cross section 
(7,?^'^^ defined through 



D 



dxpdxdQ"^ xQ'^ 



i-y+^-^]a^^'H^jp,^,Q') (7) 



D(3) 



and is related to the diffractive structure functions Fr 



and F, 



D(3) 



by 



y 



D(3) _ pD{3) _ 

" ' 1 + (1 - y f 



D{3} 



(8) 



y is defined as before, and F^'"'^^ is the longitudinal diffrac- 
tive structure function. 

The proof of Collins "7, that QCD factorization is ap- 
plicable to diffractive DIS, has the consequence that the 
DIS cross section for 7*p — > XY can be written as a con- 
volution of a partonic cross section cr^ , which is calculable 
as an expansion in the strong coupling constant as, with 
diffractive PDFs /^P yielding the probability distribution 
for a parton a in the proton under the constraint that the 
proton undergoes a scattering with a particular value for 
the squared momentum transfer t and xp . Then the cross 
section for ^*p XY is 



dxpdt ^ 



d^a2*{x,Q^0fa{^,Q^;^JP,t). (9) 



This formula is valid for sufficiently large and fixed xp 
and t. The parton cross sections are the same as those for 
inclusive DIS. The diffractive PDFs are non-perturbative 
objects. Only their evolution can be predicted with the 
well known DGLAP evolution equations, which we shall 
use in LO and NLO. 

Usually for f^{x, Q^; xp, t) an additional assumption 
is made, namely that it can be written as a product of two 
factors, fp/p{xp,t) and fa/p{P,Q'^), 

fai^^ Q^;Xp, t) = fp/p{xp, t)fa/p{P = x/xp, Q2). 

(10) 

fp/p{xp,t) is the pomeron flux factor. It gives the 
probability that a pomeron with variables xp and t cou- 
ples to the proton. Its shape is controlled by Regge asymp- 
totics and is in principle measurable by soft processes un- 
der the condition that they can be fully described by single 
pomeron exchange. This Regge factorization formula, first 
introduced by Ingelman and Schlein jl5| . represents the re- 
solved pomeron model, in which the diffractive exchange, 
i.e. the pomeron, can be considered as a quasi-real parti- 
cle with a partonic structure given by PDFs fa/p{P, Q^)- 
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P is the longitudinal momentum fraction of the pomeron 
carried by the emitted parton a in the pomeron. The im- 
portant point is that the dependence of on the four 
variables x, Q^, xp and t factorizes in two functions fp/p 
and fa/p-, which each depend only on two variables. 

Since the value of t could not be fixed in the diffrac- 
tive DIS measurements, it has been integrated over with 
t varying in the region tcut < t < tmin- Therefore we have 
according to 



f{xp) 



dtfp/p{xp,t), 



(11) 



where tcut = — 1 GeV^ and ^min is the minimum kinemat- 
ically allowed value of In the pomeron flux factor is 
assumed to have the following form 



f]p/pixip,t) — 



_ l-2ajp(t) 



exp(Bjpt). 



(12) 



ap{t) is the pomeron trajectory, ap(t) = ap{0) + a'pt, 
assumed to be linear in t. The values of Bp, ap{Q) and a'p 
are taken from ^ and have the values Bp — 4.6 GeV^^, 
ajp(O) = 1.17, and a'jp =0.26 GeV^ Usually /jp/p (a; jp, 
as written in Eq. 1)12(1 has in addition to the dependence 
on xjp and t a normalization factor A^, which can be in- 
ferred from the asymptotic behavior of crtot for pp and pp 
scattering. Since it is unclear whether these soft diffrac- 
tive cross sections are dominated by a single pomeron ex- 
change, it is better to include N into the pomeron PDFs 
fa/p and fix it from the diffractive DIS data 1 . The 

diffractive DIS cross section cr^^^'' is measured in the kine- 
matic range 6.5 < < 120 GeV^ 0.01 < /3 < 0.9 and 
10-4 <xp < 0.05. 

The pomeron couples to quarks in terms of a light fla- 
vor singlet E{zp) = u[zp) + d[zp) + s{zp) + u{zp) + 
d{zp) + s(zp) and to gluons in terms of g{zp), which are 
parameterized at the starting scale Qo = VS GeV. zjp is 
the momentum fraction entering the hard subprocess, so 
that for the LO process zp — (3, and in NLO f3 < zp < 1. 
These PDFs of the pomeron are parameterized by a par- 
ticular form in terms of Chebychev polynomials as given in 
PP . Charm quarks couple differently from the light quarks 
by including the finite charm mass rric = 1.5 GeV in the 
massive charm scheme and describing the coupling to pho- 
tons via the photon-gluon fusion process. For the pomeron 
PDFs, we used a two-dimensional fit in the variables zp 
and and then inserted the interpolated result in the 
cross section formula. 



E 



dt 



dxp / dzp 

Jo Jy^i 



dy I dx^ 



f-l/e{y)fa/-i{x^,M'^)fp/p{xp,t)fb/p{zp, Ml 

dCT(")(a6-> jets). 



(13) 



y, x~f and zp denote the longitudinal momentum frac- 
tions of the photon in the electron, the parton a in the 
photon, and the parton b in the pomeron. and Mp 
are the factorization scales at the respective vertices, and 
dcr'"' {ab jets) is the cross section for the production of 
an n-parton final state from two initial partons a and b. 
It is calculated in LO and NLO, as are the PDFs of the 
photon and the pomeron. 

The function f-y/dy), which describes the virtual pho- 
ton spectrum, is assumed to be given by the well-known 
Weizsacker-Williams approximation. 



a 
2^ 



1 + (1 
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max 



(i-y) 



2mly{ 



y 

- y 



9 2 



? 9 



)2 

- max 



(14) 



Usually, only the dominant leading logarithmic contribu- 
tion is considered. We have added the second non-logarith- 
mic term as evaluated in JHI- Qmax = 0.01 GeV^ for the 
cross sections calculated in this work. 

The formula for the cross section da^ can be used for 
the resolved as well as for the direct process. For the latter, 
the parton a is the photon and f^/^(x-^, M^) — (5(1 — x^), 
which does not depend on M^. As is well known, the dis- 
tinction between direct and resolved photon processes is 
meaningful only in LO of perturbation theory. In NLO, 
coUinear singularities arise from the photon initial state, 
that must be absorbed into the photon PDFs and produce 
a factorization scheme dependence as in the proton and 
pomeron cases. The separation between the direct and re- 
solved processes is an artifact of finite order perturbation 
theory and depends in NLO on the factorization scheme 
and scale M^. The sum of both parts is the only physically 
relevant quantity, which is approximately independent of 
the factorization scale M^ due to the compensation of the 
scale dependence between the NLO direct and the LO re- 
solved contribution [ITlTT] . 

For the resolved process, PDFs of the photon are need- 
ed, for which we choose the LO and NLO versions of GRV 
[TF| . They have been found to give a very good description 
of the cross sections for photoproduction of inclusive one- 
and two-jet final states ^9 ,10 . 



2.3 Cross Section Formula 

Under the assumption that the cross section can be cal- 
culated from the well known formulae for jet production 
in low ep collisions, the cross section for the reaction 
e+p e + 2 jets + X' + Y is computed from the following 
basic formula: 

da^{ep e + 2 jets + X' + Y) = 



3 Results 

In this Section, we present the comparison of the theo- 
retical predictions in LO and NLO with the experimen- 
tal data from HI '12'. In this paper, preliminary data on 
cross sections differential in xi!^'^'^ and z''p^ for the diffrac- 
tive production of two jets in the kinematic regions spec- 
ified in Tab. m are given. These two cross sections are the 
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only differential cross sections, which are not normalized 
to unity in the measured kinematic range. All other differ- 
ential cross sections, namely those differential in the vari- 
ables \og,,xjp, y, E^'\ Mf\ M^f , rf^'^ and \Arj^-'% 
are normalized cross sections. With these latter distribu- 
tions, only the shape can be used to test a possible fac- 
torization breaking in the resolved component. 

Before we confronted the calculated cross sections with 
the experimental data, we have corrected them for hadroni- 
zation effects. The calculated cross sections are the cross 
sections for the production of QCD jets, which consist 
either of one parton or a recombination of two partons 
according to the fcT-cluster algorithm. The experimental 
cross sections are measured with hadron jets constructed 
with the same jet algorithm. Although the difference be- 
tween the two kinds of jets is not large, in particular for 
jets with sufficiently large Et's, we have corrected the 
originally calculated cross sections with a factor Chad for 
the transformation from QCD jets to hadron jets. The 
correction factors Chad for the differential cross sections 
in the kinematic variables of interest are shown in Fig. |2| 
Here, Chad is the ratio of the respective cross sections for 
hadronic jets to partonic jets. The correction factors have 
been calculated from Monte Carlo models including LO 
cross sections together with parton showering and sub- 
sequent hadronization by the HI group As seen in 
Fig. El Chad is approximately equal to one with devia- 
tions less than 20%. The only exception is Chad for the 
^jets f,]-Qgg section with a value that is appreciably differ- 
ent from one for x^°'^^ > 0.6. 

The differential cross sections have been calculated in 
LO and NLO with varying scales, where the renormaliza- 
tion scale and both factorization scales are set equal and 
are = ^Ei^^^ with ^ varied in the range 0.5 < ^ < 2. This 
way we hope to have a reasonable estimate of the error for 
the theoretical cross sections and are not in danger to base 
our conclusions concerning factorization breaking only on 
one particular scale choice. Note that for the pomeron 
PDFs the variation of the factorization scale is restricted 
by their parameterization to Mp < 150 GeV^. 

The theoretical cross sections are presented in two ver- 
sions in LO and NLO, respectively. In the first version no 
suppression factor R is applied. It corresponds to the LO 
or NLO prediction with no factorization breaking, labeled 
i? = 1 in the figures. The second version is with a sup- 
pression factor R = 0.34 in the resolved cross section, 
labelled R = 0.34 in the figures. This particular value for 
R is motivated by the recent work of Kaidalov et al. j2Qj . 
These authors studied the ratio of diffractive to inclusive 
dijet photoproduction in the HERA regime with and with- 
out including unitarity effects, which are responsible for 
factorization breaking, as a function of Xj. In this study 
they applied a very simplified dijet production model for 
this ratio, which is very similar to the model proposed 
by the CDF collaboration for pp collisions 3- . From the 
calculations of this ratio, with and without unitarity cor- 
rections, they obtained the suppression factor R = 0.34 
for Xj < 0.3 (see Fig. 6 in Ref. [201)) which they attribute 
to the resolved part of the photoproduction cross section. 



We shall use this value of the suppression factor as a first 
try and apply it to the total resolved part in the LO cal- 
culation and to its NLO correction. The direct part is, in 
both cases, left unsuppressed {R = 1). It is clear that not 
all of the distributions will be sensitive to the value of R. 
Furthermore, most of the distributions are normalized to 
one, so that the absolute magnitude can not be used as a 
discriminator for the occurence of a suppression factor. 

Our LO (top) and NLO (bottom) results are shown 
in Fig. O for the differential cross sections in x'^*"^ (left) 

and z^p^ (right), which are not normalized to one. The 
normalized distributions in x'^*"^, z^p^ , log;^Qa;_p, y, 

M]^'", Mj2*', and |Z\?7j'=*"| are shown in LO and NLO 
in Figs.l#l 

For da/dx-!^*'^ (Fig- El left), we have very different cross 
sections for R = 1 and R — 0.34 and for the scale choice 
^ = 1. An exception is the highest x^°^^-hm, where the 
difference is only 20%, since in this bin the direct contri- 
bution is dominant and the suppression factor is therefore 
less effective. In all the other bins, da/dxi^^^ with R = 0.34 
is reduced by this factor as expected. Except for the high- 
est a;iy''*''-bin, neither of the two LO calculations agrees 
with the data. The R = 1 cross section is too large and 
the R = 0.34 cross section is too small. Only when we 
consider the scale variation with 0.5 < ^ < 2 as a realistic 
error estimate, we would conclude that the unsuppressed 
LO cross section {R — 1) is marginally consistent with 
the HI data inside the experimental errors. At NLO, the 
conclusion is reversed: the suppressed cross section now 
agrees very well with the data, while the unsuppressed 
cross section drastically overestimates the data. 

For da/dz-'p'^ in Fig. |3| (right), the agreement of un- 
suppressed and suppressed cross sections with the data is 
equally marginal at LO, even within the respective error 
bands, while it is excellent for the suppressed NLO cross 
section. We remark that the suppressed and unsuppressed 
cross sections with ^ = 1 differ approximately only by 
a factor 0.5, since in this distribution the direct and re- 
solved contributione are superimposed differently than in 

For the normalized xi^^^ distributions in Fig. 01 (left), 
the overall agreement is, of course, better. In particular, 
the unsuppressed LO distribution agrees now with the 
data within the scale uncertainty, whereas at NLO it is 
again the suppressed distribution that describes the data 
best. Furthermore, the scale uncertainty is substantially 
reduced in the normalized distributions as expected. For 
the z^p^ distributions in Fig. ^ (right), both the unsup- 
pressed and suppressed LO distributions agree with the 
data within errors, while at NLO agreement is only found 
for the latter. 

The comparison of the normalized distributions in log^Q 
X p and y is shown in Fig. O Here the theoretical predic- 
tions for R = 0.34 and R = \ differ very little. This is 
understandable, since the xp and y dependence of the 
cross section factorize (see Eq. (tH^ 'l to a large extent. 
Only through the correlations due to the kinematical con- 
straints we observe small differences between the R = 0.34 
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Fig. 2. Ratios of hadronic to partonic dijet cross sections for two different fits of the parton densities in the pomeron |1I4| as 
well as their average 



and the R ^ I cross sections, particularly in the y dis- 
tribution. From this comparison no definite conclusions 
concerning the suppression can be drawn. All theoretical 
predictions agree more or less with the data. In the high- 
est logj^QXj^ bin the measured point lies higher than the 
theoretical points. This can be explained, at least partly, 
by an additional sub-leading Reggeon contribution, which 
has not been taken into account in the diffractive PDFs 
we are using (see Fig. 7 in [T^'l. 



Next we look at the Ei^*^ distribution in Fig. |S1 The 
LO (left) and NLO (right) distributions with R = 0.34 are 
flatter than the unsuppressed distribution as we expect it, 
since the resolved component occurs dominantly at the 
smaller iJ;^*^. The suppressed cross section agrees better 
with the data points, even if the scale uncertainty is taken 
into account. Due to the normalization of the cross section, 
the differences between LO and NLO are almost invisible. 



The distributions 1/a da/dM^-^*"^ and l/a da/dM{f'' 

/ V ^''p'^^-y^'^ ■ Although 



are correlated due to Afj^ 



M. Klasen, G. Kramer: Factorization Breaking in Diffractive Photoproduction of Dijets 



7 





900 




800 




700 




600 






£? 


500 


X 




■a 


400 


"B 


■a 






300 




200 




100 



0, 



ep e+2jets+X'+Y 

HI Preliminary 
LO, R=l 
LO. R=0.34 



0.1 0.2 0.3 0.4 0.5 . 0.6 0.7 0.8 0.9 1 

jets 



900 

800 

700 

600 

„^ 500 
S, p. 
n" 

5 400 

D 
■a 

300 
200 
100 
0, 



ep e+2jets+X'+Y 

HI Preliminary 
LO, R=l 
LO, R=0.34 



0.1 0.2 0.3 0.4 0.5 .0.6 0.7 0.8 0.9 1 

jets 
Zjp 



1400 



1200 



1000 



a. 800 



X 

1 600 
■a 



400 



200 



ep e+2jets+X'+Y 



0, 



HI Preliminary 
NLO, R=l 
NLO, R=0.34 




0.1 0.2 0.3 0.4 0.5 . 0.6 0.7 0.8 0.9 1 

jets 



1400 
1200 
1000 



% 800 



ep e+2jets+X'+Y 



N 

■a 
■a 



600 
400 
200 
0, 



HI Preliminary 
NLO, R=l 
NLO, R=0.34 



0.1 0.2 0.3 0.4 0.5 .0.6 0.7 0.8 0.9 1 

jets 

Z,r, 



Fig. 3. LO (upper) and NLO (lower) cross sections for diffractive dijet photoproduction as functions of x^^^^ (left) and 
(right), compared to prehminary HI data. The shaded areas indicate a variation of scales by a factor of two around Ei^^^ . 



the distributions in a;i^'^ and z^p^ are bound to reveal more 
detailed information on possible factorization breaking, 
we have calculated the mass distributions nevertheless. 
The results and the comparisons with the data are shown 
in Fig. [71 The experimental cross sections increase with 
Afj^***, while they decrease with increasing Mi^^ This is 
due to the correlation mentioned above. The distribution 
in MI°2^ is also correlated with the distribution in Ei^^^. 
For the mass distribution of the dijet final state, which can 
directly be measured experimentally, the LO and NLO, 



suppressed and unsuppressed distributions are very sim- 
ilar and agree with the data. In contrast, the hadronic 
mass M^*'' has to be reconstructed and is very sensitive 
to systematic errors in the measured variables. The theo- 
retical prediction follows the increase in the data only in 
LO, while at NLO the dependency is reversed and is very 
sensitive to the presence of a possible third parton in the 
final state X. 

The distributions in^^"'^ and |Z\r7j'^*''| presented in Fig.|Sl 
involve a delicate superposition of direct and resolved con- 
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Fig. 4. Normalized xi'j'^^ (left) and z'p^ (right) distributions in LO (top) and NLO (bottom), compared to preliminary HI data. 



tributions. In particular, the direct (resolved) process dom- 
inates for negative (positive) rj. While the LO 77-distri- 
bution agrees better with the data, if the resolved process 
is not suppressed (i? = 1), the conclusion is again reversed 
at NLO, as was already the case for the xjy''*''-distribution 
in Fig. 121 For the lowest bin in rj, we observe an excess 
of the theoretical prediction over the data, which is well 
known from studies of inclusive jet production at the very 
low transverse momenta studied here and which can be re- 
lated to additional hadronization effects. The distribution 
in |Z\77-'''''*| is intimately linked to the angular distribution 



of the partonic scattering matrix elements. It is thus less 
sensitive to the superposition of direct and resolved pho- 
ton contributions, and the theoretical predictions agree 
almost equally well. 



In summary, we conclude that for most LO distri- 
butions the unsuppressed theory, i.e. with no factoriza- 
tion breaking, agrees better with the experimental data. 
This conclusion is, however, premature, since at NLO it 
is the suppressed theory, i.e. with factorization breaking 
and R — 0.34, which is preferred. 
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Fig. 5. Normalized logi^xip (left) and y (right) distributions in LO (top) and NLO (bottom), compared to preliminary HI 
data. 



In [201, the suppression factor of i? = 0.34 was deduced 
from a calculation of the ratio of diffractive and inclusive 
dijet photoproduction at HERA as a function of Xj for 
two cases: (i) no absorption and (ii) absorption included. 
The calculation of this ratio for the two cases was based 
on a very simplified model, in which the ratio depended 
only on the gluon PDFs of the pomeron and proton in the 
numerator and denominator, respectively. It is of interest 
to see how this ratio behaves as a function of x'^^^ for 
the two cases i? = 1 and R = 0.34 in LO and NLO in 
the more detailed theory presented in this work, i.e. in 



a theory where this ratio is calculated from the full cross 
section formula in Eq. (|13|l and the corresponding formula 
for the inclusive dijet cross section with quarks and gluons 
and realistic experimental cuts. 

The result is shown in Fig.EKleft), where we have used 
the CTEQ5M1 parameterization for the proton PDFs 
in the inclusive cross section results. In LO and for R = 1, 
the ratio da'^'^' / da"'''^ starts at small xiy"*" = 0.05 at a 
very low value (~ 0.001) and then rises monotonically 
up to 0.032 and 0.037 at cci^*" = 0.85 and 0.95. With 
R = 0.34, i.e. with suppression of the resolved part, the 
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Fig. 6. Normalized distribution in LO (left) and NLO (right), compared to preliminary HI data. 



increase of this ratio is very much reduced. It goes up to 
0.011 at xi;'" = 0.85. At xi;*" = 0.95 the ratio is substan- 
tially larger, since in this region the unsupprcssed direct 
cross section dominates. We see that up to x^J^^^ = 0.85 
the suppressed ratio {R = 0.34) is reduced approximately 
by a factor of three as compared to the unsuppressed ra- 
tio {R = 1) as expected. The behavior of the ratio is 
somewhat different for the NLO case. In particular, the 
diffractive NLO resolved contribution has a steeper rise 
at xl^^^ = 0.25 and flatter behavior above, which is re- 
flected in both the unsuppressed and the suppressed sum. 
Compared to the corresponding curves for da^^^'^ / da"^^^ 
in pn], the qualitative behavior of our curves, in LO and 
NLO, is similar. The 'no absorption/absorption included' 
curves in [211] resemble more our LO than our NLO re- 
sults as expected. We have to keep in mind, however, that 
the kinematic constraints applied in differ from ours, 
which are the same as in the experimental analysis. This 
translates mainly into a different (smaller) normalization 
of our results. Clearly it would be interesting to measure 
^^diffr^^p.inci ^ function of x^^^^ in order to have another 
observable for measuring the suppression as a function of 
x^^^^. Compared to the cross section da/dxi^^^ considered 
earlier, this ratio has the advantage to depend less on the 
photon PDFs, which appear both in the numerator and 
the denominator and should cancel to a large extent. 

It may well be that our procedure to describe the fac- 
torization breaking by applying a suppression factor to 
the total resolved cross section is not correct and must 
be modified. An indication for this is the fact that the 
separation between the direct and the resolved process 
is not physical. It depends in NLO on the factorization 
scheme and scale M^, as already mentioned earlier. The 
sum of both cross sections is the only physically relevant 



cross section, which is approximately independent of the 
factorization scale M-y. By multiplying the resolved part 
with the suppression factor R = 0.34 the correlation of the 
M-y-dependence between the direct and the resolved part 
is changed and the sum of both parts has a much stronger 
Mj dependence than for the unsuppressed case {R = 1). 
This is shown in Fig. El (right). We see the compensation of 
the M^-dependence between the NLO direct cross section 
(dotted line) and the LO resolved cross section (dashed 
line) in the unsuppressed {R = 1) case, leading to a fairly 
My independent sum of both contributions (full line) [7| 
II 7| . When the LO resolved part is suppressed with the 
factor R = 0.34, the compensation is reduced, and the 
sum of the NLO direct and LO resolved parts becomes 
much more A/^-dependent than before (although not too 
much in the range 0.5 < M^/iJr.max < 2, as seen by the 
dashed-dotted curve in the right part of Fig.|51). 

The compensation of the Af^-dependence between the 
NLO direct and LO resolved cross section occurs via the 
anomalous or point-like part of the photon PDFs. This 
means that this part of the PDFs is closely related to the 
direct cross section. It is usually assumed that the direct 
part obeys factorization and has no suppression factor. So 
the point-like part in the photon PDFs should not be sup- 
pressed either, and the suppression factor should be ap- 
plied only to the hadron-like part and the gluon part of the 
photon PDFs. Since all three parts, point-like, hadron-like 
and gluon, are correlated through the evolution equations, 
it is not clear how this suggestion could be realized. Of 
course, if the point-like part is not suppressed, the problem 
of the insufhcient compensation of the scale dependence of 
the NLO direct and LO resolved part would be solved. It 
is, however, conceivable that this problem would be solved 
quite naturally if one attempts to incorporate absorptive 
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Fig. 7. Normalized AI^-^ ° (left) and Mjj ° (right) distributions in LO (top) and NLO (bottom), compared to preliminary HI 
data. 



effects into tfie NLO tlieory following, for example, the 
work of jnj. 



4 Conclusions and Outlook 

The recent measurement of diffractive dijet photoproduc- 
tion combined with the analysis of diffractive inclusive DIS 
data in terms of diffractive PDFs offers the opportunity 
to test factorization in diffractive dijet photoproduction. 
For this purpose we have calculated several cross sections 



and normalized distributions for various kinematical vari- 
ables in LO and NLO and compared them with recent 
preliminary HI measurements ^21- In LO we found that 
the measured distributions und unnormalized cross sec- 
tions agree quite well with the theoretical results if, by a 
reasonable variation of scales, a theoretical error is taken 
into account. This means that in a LO comparison there is 
no evidence for a possible factorization breaking expected 
for the resolved contribution. However, it is well known 
that for dijet photoproduction NLO corrections are very 
important for the direct and in particular for the resolved 
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contributions to the cross section. Indeed, the theoretical 
results at NLO disagree with the data for unnormahzed 
cross sections like da/dx^°^^ and da/dz^p^ . Agreement be- 
tween data and theoretical resuhs is found, however, if the 
resolved contribution is suppressed by a factor R = 0.34. 
This factor is motivated by a recent calculation of absorp- 
tive effects in diffractive dijet photoproduction (20! ■ Since 
NLO results are more trustworthy than any LO cross sec- 
tion calculations, we consider our findings a strong indica- 
tion that factorization breaking occurs in diffractive dijet 



photoproduction with a rate of suppression expected from 
theoretical models. 

It would be interesting to investigate hard diffractive 
photoproduction of other final states, for which the su- 
perposition of direct and resolved contributions is differ- 
ent. Such diffractive photoproduction reactions are, e.^., 
large-pT inclusive singlc-hadron production, heavy-flavour 
production with or without jets, and prompt photon pro- 
duction. In order to verify that factorization breaking dis- 
appears when the of the virtual photon is increased 
from small to larger values, it would be desirable to have 
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measurenients of diffractive production of the final states 
mentioned above as a function of Q^. 

Finally factorization breaking is expected not only in 
the diffractive region, xp ^ 1, but also at larger values 
of X]p where Regge exchanges other than the pomeron 
occur. For example, pion exchange is strong in all reactions 
with a leading neutron. Here, dijet photoproduction with a 
leading neutron has been studied in LO and NLO and 
compared to ZEUS experimental data j22|. This process 
could also be a candidate for factorization breaking in the 
resolved contribution. 
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